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Abstract 

The paper presents results of investigations of small bodies dynamics in a 
vicinity of giant planets. We used the most simple body model: gravitation- 
ally bounded, rotating contact binary affected by the tidal force acting from 
a planet. Spin variations of such binaries were extensively studied during 
planetary close encounters. Two main types of dynamical behaviour were 
observed: (i) huge but interim fluctuations of the angular velocity and (ii) 
permanent changes of a rotation during a close approach. The first type is 
observed mainly for fast rotators, while the second one was encountered in 
a population of slowly spinning objects with periods longer than 12 hours. 
Conclusions on usability of such crude physical body models in dynamical 
investigations and a comparison to previous results were attached. The re- 
sults allow us to formulate a thesis explaining the phenomenon of creation 
of the extremely slow rotators and an observational excess of such type of 
objects. 
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1. Introduction 

The phenomenon of tidal interactions was widely discussed in different 
models over last 160 years. Last decade of the XX century brought important 
papers which extensively explored the problem of tidal effects in a vicinity 



Email address: r.gabryszewski@cbk.waw.pl (Ryszard Gabryszewski) 



Preprint submitted to Planetary & Space Science 



September 8, 2011 



of planets ( IChauvineau fc Farinellal . Il995l . iBottke et al.l . Il999l ). Studies of 



small bodies dynamics in such region are usually made using the most com- 
plex physical body models which are able to give output in reasonable time 
scales. Test bodies are most often conglomerates consisting of hundreds or 
even t housands small and equal spheres interacting grav itationally with each 
other ( ISharma et al.l . l2006l . IWalsh &: Richardsonl . |2006| ). also often affected 
by the forces of viscosity. Such multiparameter models describe alterations 
of the physical structure very well but studies of dynamical processes dur- 
ing planetary approaches do not always need complex models. Reduction of 
model parameters is perceived as a loss of a part of information but crude 
models also have some advantages. They allow to get general information 
on physical phenomenons much faster than complex ones and permit to find 
dependencies between model parameters and type of a body's behaviour with 
much smaller numerical effort. Another reason is the ease of implementation 
of features which are hard to apply in complex structures, e.g. an inho- 
mogeneity defined as difference in densities and/or diameters of the body's 
components. 

The main aim of this paper is to investigate dynamics of a contact binary 
during close approaches to a giant planet. A contact binary can be treated 
as a toy model of a rubble-pile body and easily allows to include an inho- 
mogeneity of its components. In contrary to recent papers we investigated 
tidal dynamics of small bodies with diameters similar to size of Centaurs and 
cometary densities. We also tried to find out how simplification of a physical 
body model affected information on dynamical evolution in a reference to 
previous works. 



2. Model and calculations 

The model of a small body consists of two spherical components held to- 
gether by autogravity forces. Both fragments contact each other and revolve 
their common mass centre. This binary system approaches a planet at a 
small but non-collisional minimum distance. Dynamics of the binary is con- 
sidered in a N-body problem with Sun and a giant planet on an elliptic orbit. 
All the bodies in the system move in only one plane - the orbital motion of 
the binary can be prograde and retrograde while the motion of a planet can 
only be prograde. The binary rotates in the plane of orbital motion of the 
planet. This means that the torque direction and angular momentum vector 
are always perpendicular to that plane. 



2 



y 




Figure 1: Model of a test body consists of masses mi and m2 approaching a planet 
Mpi. Orbital motion of all the bodies is considered in the mass centre reference frame. 
Rotation of the binary is investigated in the frame related to its mass centre: O'. Vector d 
determines distance and direction between Mpi and the mass centre of the binary. Phase 
angle (j) defines position of the binary in non- rotating Cartesian reference frame O'x'y' 
while a is an orientation angle of the binary in respect to d vector. 
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2.1. Orbital and rotational motion 

Two types of reference systems were used: orbital motion was considered 
in the mass centre reference frame of all the bodies while rotation of the 
binary was investigated in a frame related to its own mass centre, see Fig. 
[U The equation of orbital motion can be expressed as follows: 



N 



(1) 



where i=l,..,N and N is the number of bodies in the system, fij = fi — 
^j) ^ 7^ ^ is the gaussian gravitational constant, nij is the mass of jth 
body, and and fj are position vectors of ^th and jih. body respectively. 
The equ ation was solved numerically using the re current power series (RPS) 
method (IHadjifotinou fc Gousidou-Koutital . Il998l ). We extended the original 
2-body algorithm to solve N-body problem. The auxiliary variables were 
defined as follows: 



Pi ^ij Sij ^ijy ^ij ^ij^ij^ Qij 



Ui 



(2) 



where fj is the vector of velocity of the zth body. If we put them to the 
equation [H we get: 
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Pi = -k'^J2'^jij 



(3) 



If we designate the derivatives: 



2rjj ■ Tjj 
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we can determine the recurrent dependencies (written in the order of execu- 
tion) : 
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where n is the coefficient number of the recurrent power series. The procedure 
of integration wa s enriched by t he function for an automatic optimum time 
step calculation ( iSitarskil . Il979l ). The value of the time step is determined 
using last two coefficients of the recurrent power series: 



^n-l — l^^n-ll + IZ/n-ll + |^n-l| 

An = \Xn\ + \yn\ + \Zn\ (6) 

where x, y and z are rectangular coordinates of a body (z coordinate is zero 
in the case considered in this paper). Then time step can be defined as: 



^loc 

An-1 + Anhoil + 



(7) 



where: 



and Eioc is the wanted accuracy of numerical solution. The formula ([7]) allows 
us to determine the time step for every body in the integration process, we 
should choose the smallest value. Positions and velocities of the zth body 
can be obtained from the equations below: 
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r1 = j:^jh^-\ (8) 

n 

i=i 

Variations of a binary's angular velocity were investigated using Euler 
equations. In the planar problem where we have only one component of the 
angular velocity, we can simplify these equations to the form: 

D = luj (9) 

where / is the total moment of inertia of a binary, D is its torque and u is 
the angular acceleration of the body. The moment of inertia is defined as a 
sum of moments of both fragments: 

I = j2{lmid^ + mirf) (10) 
1=1 ^ 

where mi are the masses of conglomerate's components, di are the radii of 
mi, and ri are the lever arms. The total torque is defined as follows: 

2 

D = Y,fiXFMm, (11) 
1=1 

where FMmi are forces acting from a planet M to masses mi respectively. 

The equation was solved numerically by the Runge-Kutta method of 
the 5th order (RK5) with automatic time step calculations using standard 



algorithms published in the Numerical Recipes (iPress et al.l . Il992l ). Minor 
modifications were made to get a double precision output. Two different 
time scales of orbital and rotational motion were connected by invoking the 
RK5 integrator on every substep of the RPS method. 

2.2. Breakup of a contact binary and reconnection of the fragments 

We defined the tidal disruption as a separation of fragments due to tidal 
force. This separation is accelerated by a centrifugal force caused by the 
binary's rotation. Fragments start to be treated as separate bodies when the 
difference of accelerations acting from the planet on binary's components - 



6 



increased by the centrifugal acceleration - exceeded the mutual acceleration 
in the mim2 system: 



a < \ai — 02! + flc/ (12) 

where a is a gravitational acceleration of nii respectively to m2, ac/ is an 
acceleration due to centrifugal force and ai and a2 are accelerations of mi 
and 777-2 respectively caused by gravitational force acting from the planet. 

When fragments are separated their angular momentum do not change in 
time because of their spherical shape. The fragments can be merged again by 
mutual collisions. We assumed collisions to be pure inelastic in our model, 
this means that every mutual collision of fragments creates back a contact 
binary. This body can be disrupted in the next step of an integration process 
if the condition f lT^ is still fulfilled. These assumptions allowed to calculate 
new angular momentum and angular velocity of the body in moments of its 
creation and disruption. 

2.3. Initial conditions 

The investigated systems consisted initially of 3 bodies: Sun, a planet 
(Jupiter and Neptune in the calculations) and a contact binary. We gen- 
erated 1400 sets of initial conditions. Starting perihelia of contact binaries 
were randomly spread in a distance to 0.003 AU from the perihelion of a 
planet, and the aphelia were defined between 9 and 45 AU: the bodies (a 
planet and a contact binary) started the dynamical evolution close to their 
perihelia. Their mean anomalies had initial values defined between 357 and 3 
degrees (which defined initial distances between planet and the centre of the 
mass of the binary between 0.1 and 0.9 AU). Values of orbital inclinations 
of binaries were set to and 180 degrees - prograde and retrograde motion 
was investigated. Masses of contact binaries were parametrized by densities 
and diameters of the fragments. The ranges of fragments diameters were 
defined between 50 and 100 km which interacts ver y well with diameters o f 



the largest Centaurs: 2060 Chiron and 5145 Pholus ( IStansberry et al.l . l2008l ) 



The ranges of de nsities were chosen very wide, fro m low comet ary 0.2 g /cm' 
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( ISekaninal . 119821 ) to Centaur values of 1.5 g/cur ( IGrundya et al.l . 120071 ). 

Most of previous papers studied the problem of tidal encounter in a 2- 
body problem where the small body moved on parabolic or hyperbolic orbit 
with a planet in its focus. Initial conditions were determined by small body's 
velocity "at infinity" Vqo and by periapse distance q. In our initial conditions 
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varied in a range between 3 and 4 km/s for prograde motion and 30 to 34 
km/s for retrograde motion. The minimum distance between a planet and a 
binary could only be known a'posteriori. 

We investigated variations of orbital elements of all the bodies in the 
system, angular velocity and geometry of the contact binary together with 
variations of mutual distance between its fragments after the breakup. Equa- 
tions of motion were integrated up to 6 subsequent planetary encounters. 

3. Results 

The simulations allowed to distinguish two frequently observed types of 
behaviour: oscillations around initial value in case of fast (above 4 rev/day) 
and moderate (~2.5 rev/day) rotators and permanent changes of the angular 
velocity of slow rotators (see Fig. [2]). Small bodies with starting angular 
velocity of 0.1 rev/day or lower can even change their spin directions. The 
centrifugal force can be an important factor which helps tides to disrupt 
fast rotators, see the evolution of bodies with initial spin of 3.5 and 4.5 
rev/day (the self disruption due to sole centrifugal force is about 5.2 rev/day 
in presented case). The reconnection of fragments is stable if the spin velocity 
of newly created binary is much lower than the initial one (see the case 
of a;=4.5 rev/day). Distinct variations of a binary's spin were observed in 
distances up to 10 planetary radii for every investigated giant planet. We 
noticed that the influence of tidal force was nearly insignificant in a distance 
above 25 planetary radii (standard values of planetary equatorial radii were 



The difference in spin variations between fast and slow rotators are the 
derivative of the torque D changes. The total torque changes its sign 4 
times while an aspherical body revolves around its spin axis. This means the 
resultant torque is usually small for fast rotating objects but for slow rotators 
the change of torque is able to achieve much larger values. Modifications of 
angular velocities are done in a quite short time: 12 to 24 hours. When a 
distance to a planet is greater than 0.01 AU, torque makes only insignificant 
modifications to angular velocity - changes seem to be random and do not 
accumulate during subsequent close-ups. 

The contact binary body model allowed to incorporate an inhomogeneity 
in a very simple way. We made lots of simulations with fragments of different 
mass ratio (assuming the same sum of the masses and fragments densities). 
The general observed dependence between evolution of angular velocity and 
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Figure 2: Variations of angular velocity of a contact binary in a function of time for 

different starting values. Upper figure shows the distance between Jupiter and the small 
body. Initial distance between the planet and the body is 0.197 AU, while the minimum 
distance of the approach is 0.0034 AU (~7 planetary radii). Solid flat line on the upper 
chart denotes the planetary radius. The curves are interrupted when the binary breaks up. 
After the break-up event fragments are able to adjoin (case of 01=4.5 rev/day) creating a 
new contact binary with much lower angular velocity. 
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mass ratio of fragments is: the larger disproportion between masses, the 
lower amplitude of velocity variations, see Fig. [31 Distinct difference in the 
amplitude of spin velocities was observed if masses of fragments differed 
more than 15 % . Such dependence is an effect of the body's shape. If the 
fragments have comparable masses and sizes, the shape of the body is the 
most elongated, contact binary formed by such fragments can be treated as 
an approximation of the most possible aspherical body created by 2 spheres. 
The influence of torque D is the largest in such and gives us in effect 

largest possible variations of spin velocity and distance of disruption. 

Fig. m shows the maximum distances of tidal disruption for contact bina- 
ries with densities: 0.2 g/cm^ (lower graph) and 0.6 g/cm^ (upper graph). 
Both graphs include curves showing the maximum distance of disruption for 
non-rotating binary (where tidal interactions only were taken into account) 
and rotating with spin velocity of 2.67 rev/day (tidal interactions + rotation 
P=9h). The figure shows that distance of disruption depends mainly on frag- 
ments mass ratio and their densities. Fast rotation can significantly increase 
this distance only when the binary has a very low density and consists of 
fragments of similar masses. 

Fig. [5] shows variations of angular velocity in a function of a phase angle 
0. The computations were made for different starting values of angular ve- 
locity but the figure shows only the cases oi u = 1.5 and 3.5 rev/day. Two 
main types of the behaviour were observed. Below some threshold value 
(2.5 rev/day in the presented case) rotation is slowing down for most values 
of the (f) angle, the change of the spin direction is also possible (see lower 
graph). Although for some specific range of the angles - 100 to 150 de- 
grees - angular velocity could even be increased up to ~30 per cent. That 
geometry-dependent effect is related to an elongated shape of a contact bi- 
nary and the changes of the torque. The upper graph shows that fast rotators 
are able to increase their angular velocity regardless of a starting value of 
0. If condition (fT2l) is fulfilled, the binary disrupts and fragments start to 
evolve separately. Time of disintegration depends on the phase angle. In 
some cases fragments can adjoin and create a new binary rotating with an 
angular velocity much lower than the initial one. 



Tidal force can be responsible for observed excess of slow rotators (iFarinella et al 



198ll ) but there are als o other processes a s like non-tidal brea kups, non-tidal 



splitting flHarrisl . |2002|) and YORP effect flBottke et al.l . l2006f ) which can also 
create slow rotating bodies. Our results show that tidal interactions can be 
an effective process of slowing down low density objects in a vicinity of giant 
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Figure 3: Variations of contact binary's angular velocity in a function of time during close 
encounter to Jupiter. Upper figure shows the diagram for identical masses of fragments, 
lower figure is for mass ratio of 95/5 (assuming the same sum of masses and the same 
fragments densities in both cases). The curves are interrupted when the binary breaks up. 
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Figure 4: Maximum distance of tidal disruption in a function of fragments mass ratio 
for non-rotating contact binary (when sole tidal interactions were taken in account, tides 
only description on charts) and rotating with angular velocity of 2.67 rev/day (tides + 
rotation P=9h). Binary's densities were set to 0.6 g/cm"^ (upper chart) and 0.2 g/cm^ 
(lower chart). Solid flat line on the figures denotes the planetary radius of Jupiter. 
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Figure 5: Variations of a contact binary's rotation as a function of time. Starting value 
of angular velocity is set to 3.5 rev/day on upper graph and 1.5 rev/day on lower graph. 
Fast rotators can break up close to a minimum distance of planetary approach, the time 
of disruption depends on the phase angle 4>. After the break-up event fragments are able 
to adjoin (cases of 0=130 and 150 deg) creating a new contact binary which rotates with 
much lower angular velocity. On lower graph most of slow rotators tend to lower their 
angular velocities after the close approach. Binaries with starting (p between 10 and 50 
deg changed their spin directions permanently. 



planets - in distances close to 7 - 8 planetary radii and smaller. In contrary to 
Weidenschilling's tidal despinning (1989), tidal interactions with giant plan- 
ets can be treated as the mechanism able to produce objects with extremely 
long r otational periods, even co mparable to the one of 288 Glauke - over 50 
days (IKryszczyhska et al.l . l2003l ) . 

Fast rotators are far less susceptible for permanent changes of angular 
velocity but such type of dynamics was also observed. The increase of the 
spin velocity leads to tidal disruption in most studied cases. 

The results presented on figs. |2l |3] and [5] concern contact binaries with 
very low densities of 0.2 g/cm^. But general dependencies for the bodies with 
densities of 0.6 and 1.5 g/cm^ are very similar. The differences concern the 
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distances of disruption which are 1.5 to 2 times smaller respectively. Decrease 
of the breakup distance strongly affects variations of angular velocity - much 
higher amplitudes are observed. The centrifugal force hardly increases the 
disruption distance for bodies with densities of 0.6 g/cm^ and greater (see 
Fig. SD. 

About 33% of tested binaries were disrupted during the planetary ap- 
proaches. In 72% of those cases, one of the disrupted fragments hit the 
planet while another passed it by on a highly eccentric orbit. Both frag- 
ments survived the planetary close-ups in about 27% cases. We also found 
very peculiar dynamics: double systems were created in about 1% of cases. 
The separation of the fragments varied from tens of meters to about 180 
km. Orbits of such double system had very elliptic shape: their eccentrici- 
ties varied between 0.1 to 0.6 with higher values mainly for compact system. 
Such system were very unstable, they were disintegrated in a few to several 
days. The rest of double systems were stable in a timescale of computations 
but their mass centres moved on hyperbolic orbits which means they will 
leave the planetary system. Those peculiar double systems were created by 
collisions which were defined to be pure inelastic in our model. 

4. Conclusions 

The conclusions can be summarized as follows: 

• We can distinguish two main types of behaviour of contact binaries 
angular velocities approaching giant planets: 

(i) large but temporal amplitude of variations during a close approach 
without significant modification of the initial value after encounter, (ii) 
permanent change of a spin velocity of slow rotators with periods longer 
than 12 hours. 

• Modifications of angular velocity are made in 12 to 24 hours during the 
close up. 

• Binaries with angular velocities of 0.1 rev/day or lower can change the 
spin direction. 

• Significant variations of angular velocities can be observed up to ~10 
planetary radii due to tidal interactions, variations above this limit are 
random and do not accumulate. Tidal force at distances of 25 planetary 
radii or greater have a very weak influence on angular velocity. 
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• Distance of tidal disruption depends on fragments densities and their 
mass ratio. Centrifugal force is able to extend this distance only in case 
of very low densities. 

• Extremely long period of rotation can be an effect of deep encounter 
between an aspherical small body and a giant planet. Spin periods of 
slow rotators can even decrease to months {u < 0.01 rev/day). 

4.I. Discussion of results 

In Sharma's paper (2006), Figs. 8 and 10 show variations of angular 
velocity in a function of time during planetary close-ups. Although the data 
are in non-dimensional units we can simple convert them to our units using 
the information placed below the diagrams. The body which is initially a 
fast rotator (a;=4.4 rev/day) is able to decrease its spin velocity to a value of 
u)=1.84 rev/day due to tidal interactions. Initial moderate or slow rotators 
are able both to substantially increase their periods to a;=4.52 rev/day or 
slow it down to u=0.6 rev/day. These levels of variations are in a very 
good agreement with the results presented in this work. We confirm that 
fast rotators can also essentially decrease their spin velocity but this type of 
evolution is observed seldom in our calculations. 

The time of spin rate changes is about 7 hours at Sharma's results com- 
paring to 12 - 24 hours in our outcomes. This difference can be explained by 
the minimum planetary distance during the approach: 1.8 planetary radii at 
Sharma's paper vs 7 planetary radii in our outcomes (voo was about 3 km/s 
both in Sharma's and our calculations), and also by a difference in the mass 
of an approached planet. 

Unfo rtunately, we cannot d o similar comparison to most of previous 



pape rs (iBottke &: Meloshl . Il996l . IScheers et al.l . |2004 IWalsh fc Richardsonl . 



20061 ) . The authors described mainly the results of disruptions, i.e. the num- 
ber of created fragments and their general dynamical properties, the works 
do not include numbers for angular valocity variations during planetary ap- 
proach. But general dependencies in this work are fully consistent with the 
outcomes of the papers of Walsh & Richardson (2006), related to test bodies 
with flattened ellipsoidal shape. 

It should be noticed that Solar System model used in Sharma's and 
Walsh's papers were completely different comparing to assumptions made 
in this work: their calculations were made in a 2-body problem where small 
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body approached a planet on a hyperbolic orbit. But their dynamical con 
elusions are nearly the same a s ours. The presented outcome is also con 



sistent with Sekanina's works (ISekanina fc Yeomand . Il985l ISekanina et al. 



I994J . 119981). In those papers they considered the importance of centrifu- 
gal force as a factor responsible for disruption of low density (0.2 g/cm^) 
comets during close encounters with a planet. Their results also stay in good 
agreement with our outcomes. 

The results presented in this work show that detail investigations of ob- 
jects dynamics can be driven using quite crude physical body models if im- 
portant dynamical processes are included. 
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